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Accreting magnetars: a new type of high-mass X-ray binaries? 
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The discovery of very slow pulsations (f spin = 5560>s) has solved the long-standing ques- 
tion of the nature of the compact object in the high-mass X-ray binary 4U 2206+54 but has 
posed new ones. According to spin evolutionary models in close binary systems, such slow 
pulsations require a neutron star magnetic field strength larger that the quantum critical value 
of 4.4 X 10'^ G, suggesting the presence of a magnetar We present the first XMM-Newton 
observations of 4U 2206+54 and investigate its spin evolution. We find that the observed 
spin-down rate agrees with the magnetar scenario. We analyse ISGRl/ INTEGRAL observa- 
tions of 4U 2206+54 to search for the previously suggested cyclotron resonance scattering 
feature at ~30 keV. We do not find a clear indication of the presence of the line, although 
certain spectra display shallow dips, not always at 30 keV. The association of these dips with 
a cyclotron line is very dubious because of its apparent transient nature. We also investigate 
the energy spectrum of 4U 2206+54 in the energy range 0.3-10 keV with unprecedented detail 
and report for the first time the detection of very weak 6.5 keV fluorescence iron lines. The 
photoelectric absorption is consistent with the interstellar value, indicating very small amount 
of local matter, which would explain the weakness of the florescence lines. The lack of matter 
locally to the source may be the consequence of the relatively large orbital separation of the 
two components of the binary. The wind would be too tenuous in the vicinity of the neutron 
star. 

Key words: X-rays: binaries - stars: neutron - stars: binaries close - stars: emission line. Be 
- stars: magnetars 
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1 INTRODUCTION 

The X-ray source 4U 2206+54 is a peculiar member of the high- 
mass X-ray binary systems (HMXB). In this type of X-ray bina- 
ries, the normal stellar component is a massive star, usually an O 
or B star, while the X-ray emitting component is generally a neu- 
tron star. Mass transfer in HMXBs may occur via the strong stel- 
lar wind of the optical companion as in systems containing super- 
giant stars, the so-called supergiant X-ray binarie^ (SGXB) or via 
the disc-like envelope of photospheric material around the equa- 
tor of main-sequence companions. In this latter case, the system is 
called a Be/X-ray binary (BeXB). Wind-fed systems are normally 
persistent X-ray sources with X-ray luminosities of the order of 
10'^ - 10^* erg s"', while most BeXB are transient sources with 
episodes of violent eruption in X-rays that may reach 10^^ erg s"' . 



* E-mail: pau@physics.uoc.gr 

' Accretion via an accretion disk is believed to occur in only three SGXBs 

due to Roche lobe overflow. 



The optical/X-ray properties of 4U 2206+54 do not fit in any of 
these two categories of HMXBs. 

It is classified as a HM XB because it contains a n eutron star 
orbiting a 0-type companion JNegueruela & Reia200ll) . However, 
4U 2206+54 does not belong to the supergiant subgroup because 
the optical companion is a main-sequence star (Blay et al. 200Q). 
The absence of any obvious long-term trends in the evolution of the 
Ha line spectral parameters and the lack of any correlation between 
the Ho- line equivalent width and infrared magnitudes and colours 
differ from what it is expected in a typical BeXB. Optical and in- 
frared studies indicate that the optical companio n in 4U 2206+54 
is an 09.5V star located at a distance of 2.6 kpc JBlav et al.ll2006l) 
with a higher than normal He abundance. The X-ray emission is 
the re sult of accretion from an abnormally slow wind (IRibo et al.l 
l2006l) . 

For a long ti me the nature of the compact companion was 
debated. Recently, iReig et alj (2009|) discovered a 5560-s modu- 
lation in the X-ray light curve that was interpreted as the spin 
period of the neutron star. The discovery of very slow pulsations 
has solved the controversy on the nature of the compact object in 
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4U 2206+54, but has prompted new questions. According to spin 
evolutionary models in close binary systems, such slow pulsations 
require magnetic field strengths of the order o f ^10'"^ G, that is , 
they require the system to harb our a magnetar ([Finger et al.lbOlCI : 
llkhsanov & Beskrovnavall201(]|) . 

The magnetic field str ength can be estimated considering the 
neutron star spin evolution. iFinger et al.l (12010) derived an average 
rate of the spin frequency change of v = (-1.7 ± 0.3) x lO"'"* Hz 
s"' . Based on the slow pulsation period and the spin down rate, they 
show that the 5560 s pulsations in the X-ray flux of 4U 2206+54 
can be explained "provided the neutron star in this system is a mag- 
netar whose surface field at the present epoch exceeds 10'"* G". 

However, these results are at odds with the tentative detec- 
tion of a cyclotron resonance scattering feature (CRSF). Although 
the statistical significance of this feature is only marginal in any 
single one measure ment, it has be en reported by thre e diflFerent ob- 
serva tories: RXTE aTorreion et alll2004.). BeppoSAX llMasetti et al.] 
l2004h . and INTEGRAL Jb lav etal J l2()o3 ; IWand 120091) at diflFer- 
ent epochs. However, no sign of this feature has been detected 
in other observat ions with the sam e instrumen ts on board RXTE 
dReig et alJlJOOgj) and INTEGRAL ( Wandl20ld) . If real, the detec- 
tion of a CRSF at ~ 30 keV would imply a magnetic field strength 
(S/IO'^G ~ £,yci/ll-6keV) of 3.3 x 10'^ G, typical of accretion- 
powered pulsars. 

In this work, we analyse the first XMM-Newton observations 
of 4U 2206+54. We confirm the existence of the 5560-s X-ray pul- 
sations, derive a new value of the spin-down rate and perform a de- 
tailed spectral study of the soft emission with unprecedented sensi- 
tivity. We also analyse ISGRl/INTEGRAL observation to search for 
the presence of the cyclotron line. 



Table 1. XMM-Newton observation log. 



Instrument 


Mode 


Filter 


Exposure (ks) 


EPIC-PN 


large window 


medium 


74.8 


MOSl 


small window 


medium 


72.4 


M0S2 


small window 


medium 


72.4 


RGS 


standard 


- 


76.4 


OM 


image 


grisml 


71.7 



spectra and light curves were generated by excluding the core of the 
PSF. The same extraction region was used for both light curves and 
energy spectra. We extracted events from an annulus region with 
inner and outer radii 10 and 60 arcsec, respectively. Background 
files were extracted from a circular source-free region with a ra- 
dius 30 arcsec centered approximately at the same distance as the 
readout node (RAWY) as the source region. To achieve the best en- 
ergy resolution and minimise further the effect of pile up the final 
PN spectrum was created from single events only. The SAS tasks 
rmfgen and arfgen were used to generate source-specific response 
matrix files and auxiliary response files for spectral analysis. We re- 
binned the EPIC-PN spectrum to oversample the full-width at half 
maximum (FWHM) of the energy resolution by a factor 3, and to 
have a minimum of 25 counts per bin to allow the use of the x^ 
statistic. 

The SAS task rgsproc was used to produce calibrated RGS 
event lists, spectra, and response matrices. The individual spectra 
of the same spectral order obtained from each one of the two RGS 
instruments were combined (rgscombine) to produce the final RGS 
spectrum. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 XMM-Newton observations 

4U 2206+54 was observed by XMM-Newton on 2011 February 
6 during revolution 2044. The observation (ObsID 0650640101) 
started at 03:39 hr UT and lasted for ~75 ks. The observation de- 
tails for the instruments on-board XMM-N ewton are presente d in 
Table [U The XMM-Newton Observatory djansen et ai]|200ll) in- 
cludes three 1500 cm^ X-ray telescopes each with an European 
Photon Imaging Camera (EPIC) at the f ocus. Two of the EPIC 
imaging spectromet ers use MPS CCD s dTumer et al.l 1200 11) and 
one uses PN CCDs (Struderet al."200lV Reflection Grating Spec- 
trometers (RGS den Herder 2001) are located behind two of the 
telescopes while the 30-cm optical monitor (OM ) instrument has 
its own optical/UV telescope (Mason et al.ll200lh . All instruments 
were in operation during the observation. Data were reduced using 
the XMM-Newton Science Analysis System (SAS version 1 1.0). 

We first reprocessed Observation Data Files (ODFs) to ob- 
tain calibrated and concatenated event lists. The EPIC-PN analy- 
sis was performed on this ODF event file. We restricted the use- 
ful PN events to those with a pattern in the range to 4 (single 
and doubles) and complying with the more strict selection criterion 
FLAG=0, which omits parts of the detector area like border pixels 
and columns with higher offset. 

Although the observation were made in "partial window" 
mode, the brightness of the source was so high that the EPIC in- 
struments were affected by pile-up. Because the PN camera has a 
faster read-out time, it is less affected by this effect. Consequently, 
we used PN data in our timing and broad-band spectral analysis. PN 



2.2 INTEGRAL observations 

INTEGRAL (INTERnational Gamma-RAy Laboratory) is an ESA 
space mission with contributions from NASA and Russia dedi- 
cated to the observa tion of the hard X-rays/soft-y-ray Universe 
dWinkler et al.ll2003r) . The two main instruments on board INTE- 
GRAL are IBIS (Imager on Board Integral S pacecraft) which op- 
erates in the 15 keV-10 MeV energy range (lUbertini et al .1120030 
and SPI (SPectometer on boar d Integral) which operated in the 
20 keV-8 MeV energy range (Vedrenne et al. 2003) . IBIS, in turn, 
is composed of two detector layers, ISGRI (Integral Spacecraft 
Gamma-Ray Imager) operational in the 15-1000 keV range, and 
PCSiT (PixeUated CSi Telescope), working in the 0.175-10 MeV 
range JLebrun et al.ll2003ULabanti et al.ll2003l) . At lower energies, 
the payload of INTEGRAL consists of the X-ray monitor JEM -X, 
sensitive to X-rays in the 3-35 keV band dLundet al.ll2003[) and 
the optical monitor OMC that covers the visual range 5000-6000 A 
(Mas-Hesse etal. 2003). 

INTEGRAL observations normally consist of a series of 
pointings called Science Windows (SCW), 2 to 4 ks each. Since the 
launch of INTEGRAL in October 2002, 4U 2206+54 has repeat- 
edly appeared in the field of view of its instruments. However, the 
number of observations with a significant detection of the source 
(those with a detection level above 7cr) is relatively small. We have 
analysed all publicly available ISGRI observations from the begin- 
ning of the mission up to revolution 966 (September 2010). We 
selected the SCW where the source was within 10 degrees from 
the centre of the pointing and the detection level was above 7. This 
selection resulted in 37 SCWs. For each SCW an ISGRI spectrum 
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Figure 1. Upper panel: PN light curve 4U 2206+54. Time zero is MJD 
55598.17 and bin size is 100 seconds. Bottom panel: FFT using a 10-s 
binned light curve. A peak at ~ 0.17 mHz, representing the 5590-s pul- 
sations, is clearly seen. 



was obtained. These spectra were combined together to produce an 
average ISGRI spectrum. The Integral Science Data Center (ISDC) 
software OSA (Offline Scientific Analysis) version 9.0 was used in 
the data processing. 



3 THE SPIN PERIOD EVOLUTION 

Prior to the determination of the spin period, we converted the ar- 
rival times to the solar system barycenter. Then we generate a light 
curve in the 0.2-12 keV energy range with a bin size of 10 seconds. 
The last ~2300 seconds of the light curve were affected by high 
background. After removing this part the total duration of the light 
curve was 72.5 ks. 

Figure [T] shows the 0.2-12 keV PN light curve rebinned to 
100-s bins and the Fast Fourier transform (FFT) obtained from 
the 10-s binned light curve. A coherent periodicity is visible in 
the light curve and confirmed by the FFT. The maximum power 
is distributed around ~0.0001775 Hz, corresponding to a period of 
~5635 s. The power spectrum of 4U 2206+54 is characterised by 
strong red noise at low frequencies, whic h is well d escribed by a 
power-law model with a slope of 1.5-1.7 dReig et al.iii2009 ). To as- 
sess the statistical significance of this peak and give a confidence 
detection limit (dashed line in Fig. [TJ, we first corrected for the 
presence of this type of noise by fitting a power-law function to the 
continuum. The power spectrum shown in Fig. [T] was created by 
multiplying the original power spectrum by 2 and dividing by the 
power-law model. 

The small number of pulse cycles covered by the light curve 
(~13) prevented us from using the pulse phase connection tech- 
nique (see e.g. IStaubert et al.ll2009l . and references therein) to de- 
rive an accurate value of the spin period because the pulse pro- 
files would not have sufficient statistics. Nevertheless, we can use 
other less accurate techniques to refine the value of the spi n period. 
We u sed the CLEAN (Roberts et al. 1987), Lomb-Sca rgle JScargle 
1982 ) and Phase Dispersion Minimization (PDM IStellingweri 
19781) algorithms. All these algorithms are implemented in the pro- 
gram PERIOD (version 5.0), distributed with the Starlink Soft- 
ware Collection. Table |2] gives the results of the timing analy- 
sis. Another technique that achieve s smaller errors is epoch fold- 
ing JLeahv et alJll983l ; lLeahvlll987l) . We used the task ef search of 
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Figure 2. Spin-down evolution. The lines represent the best linear fit to all 
the data (solid line) and after removing the oldest data point (dashed line). 



the XRONOS package and equation (6a) in lLeahvl h987h to esti- 



mate the error. We obtained P., 



5591 ± 3 s. One more deter- 



mination of the spin period was d one with the software package 
PERIOD04 (Lenz & Bregenl2005h . which uses a discrete Fourier 
transform (DFT) algorithm and allows the estimation of uncertain- 
ties through Monte Carlo simulations. This method gave 5600± 5 s, 
where the eiTor was estimated from the distribution of frequencies 
after running 1000 simulations. The mean spin period obtained as 
the average of the six methods is 5585 ± 30 s, where here the error 
is the standard deviation of the measurements. A weighted mean 
gives 5593 ± 3 s. However, it should be noted that these calcula- 
tions apply in the case that, in the absence of pulsations, the only 
noise present in the data is that due to counting statistics, i.e, white 
noise. Be cause 4U 2206+5 4 shows stron g red n oise in its power 
spectrum C Reig et al.l 120091 ). iFinger et alj ( 120101) used a more re- 
alistic approach and assumed that the natural log of the rates is 
(as opposed to just the rates) Gaussian distributed. Then they used 
the Cash statistics and compared it to the results obtained with the 
Lomb-Scargle method. They found that the pulse periods estimated 
from the peak of the Lomb-Scargle periodogram were shifted by 1- 
20 s (depending on the initial light curve) with respect to the maxi- 
mum likelihood fit (minimisation of the Cash statistic). Taking into 
account these statistical considerations we adopted a final value of 
P = 5593± 10 s. 

Figure |2] shows t he spin-down evolut ion of 4U 2206+54. The 
data points come from lFinger et al.l ( 120101 Be ppoSAX & S uzaku), 
iReig et al.l J2009l INTEGRAL & RXTE) and lwand J20ld . INTE- 
GRAL). Using only the BeppoSAX and Suzaku measurements of 
the spin period, I Finger et al.l ( 120101) obtained a spin-down rate of 

V = (-1.7 ± 0.3) X 10"''' Hz s"'. Figure |2] shows that this trend 
continues up to present. A linear fit to all the data points gives 

V = (-1.5 ± 0.2) X lO"'"* Hz s"' with a coiTelation coefficient of 
0.97 (solid line). 



4 THE XMM-NEWTON SPECTRUM 

The XMM observations allows us to investigate the low-energy part 
of the spectrum in 4U 2206+54 with unprecedented sensitivity. 
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Table 2. Spin period detennination. 



Method 


/'spin (S) 


CLEAN 

Lomb-Scai-gle 

PDM 


5580±107 
5580±107 
5532±105 


DFT 


5600±5 


Epoch folding 


5591±3 


Mean 
Weighted mean 


5586±30 
5593±3 


Final adopted value 


5593±10 



Table 3. Best-fit spectral parameters to the EPIC-PN spectrum. The lumi- 
nosity and blackbody emitting radius were computed assuming a distance 
of 2.6 ± 0.3 kpc (Blay et al. 2006). Quoted eiTors are at 90% confidence 
level for a single parameter. ABS: photoelectric absoiption, BB: blackbody, 
BMC: Bulk Motion Comptonization, PL: power law. 



ABS*(BB+BMC) 


Nh (xlO^' cm"^) 




2.23±0.02 


iT-bb (keV) 




0.51±0.02 


i?bb (km) 




1.3±0.2 


kT,„l (keV) 




J 47+0.07 
1 07+0:05 

'■^'-0.4 


a 




Lm (X 10^5 erg s-l) 




0.13 


iflMcCx 1035 ergs-') 




2.7 


^2/dof 




0.98/197 


ABS*(BB+PL) 


Nh (xlO^' cm"^) 




3.72±0.08 


<:rbb (keV) 




1.63±0.03 


Rhh (km) 




0.37±0.04 


r 




0.94±0.03 


^^norm 




1.11 ±0.04 


Lift (X 10^5 ergs-') 




1.17 


^2/dof 




1.25/200 


ABS*(BB+BB) 


Nh (xlO^' cm-^) 




2.00±0.05 


iTbbl (keV) 




0.69±0.02 


ffbbi (km) 




1.1±0.1 


kTbbi (keV) 




2. 11 ±0.02 


i?bb2 (km) 




0.34±0.04 


x^im 




1.45/200 


io.3-2keV (xlO" ergs" 


-') 


0.16 


i-2-lOkev(xWergs- 


') 


2.21 


a:xlO--phkeV-'cm- 


^s' at 1 keV 



B 0.0001 




Energy (keV) 

Figure 3. EPIC-PN energy spectrum. The BMC (dashed line) dominates the 
description of the continuum spectrum at high energies while a blackbody 
(dotted line) has been added to connect the excess at low energies. 



when the second component is represented by a power law, by Bulk 
Motion Comptonization (BMC lshrader & TitarchukI 19981) or even 
by another blackbody. The best-fit parameters of these two combi- 
nation of models are presented in Table |3] 

The fit with the reduced x^ closest to 1 is obtained with the 
combination that includes the BMC. In this model, the power law 
spectral index a yields a value larger than 1 which means that the 
Comptonization efficiency is high. Consequently, the illumination 
factor /, defined as the number of photons multiply scattered with 
respect to the number of injected photons by the thermal com- 
ponent, is >> 1, that is to say, the spectrum appears completely 
comptonized. The size of the seed soft photons production site 
can be estimated from the luminosity of the soft photons source 
Lq = ^.^.^ where L is the resulting com ptonized luminosity and 
jco = 2.SkT^gi/kTg JTitarchuk et al.lll997l) . Using the values in Ta- 
ble [3] and assuming a comptoniz ing cloud of kT^ ~ 50 keV (as 
obtained in previous works, e.g. iTorreion et all 120041) we obtain 
Lo = 1.1 X 10^5 ergs s"'. If the soft photon source is radiating like 
a blackbody of area nR^ then we have 



i?w = 0.6-^^*(kT,aJl keWy- [km] 

which yields i?w ~ 1 km. This is consistent with the size of a hot 
spot over the neutron star surface. On the other hand, if we do not 
assume any temperature for the comptonizing cloud and use instead 
the total luminosity for the BMC component, the resulting radius is 
1.6 km, still very small and compatible with a hot spot. 

Alternatively, the size of the emitting region can be estimated 
from the blackbody component. The luminosity of the soft excess 
turns out to be Li,t, = 1.3 x lO^"* ergs s"' which, for the temperature 
given in Table [3] yields a circular radiu^ of Ri,/, « 2.4 km. This 
is again too small to be produced by any material surrounding the 



4.1 The EPIC-PN spectrum 

Figure [3] shows the EPIC-PN spectrum of 4U 22064-54, covering 
the energy range 0.2- 1 2 keV. This spectrum can be well represented 
by two-component models, affected by photoelectric absorption. 
The first component is a blackbody, which accounts for the softer 
part of the spectrum, while the second component dominates above 
~3 keV. Good fits to the EPIC spectral continuum are obtained 



^ If the interpretation that the blackbody component represents emission 
from hot spots is correct, then it seems more natural to assume a circu- 
lar area as the emitting region. However, the values found in the literature 
normally assume a spherical emitting region. For this reason, the radii in 
Table |3] correspond to a spherical emitting area. Note that if the soft pho- 
tons source is taken to radiate as a blackbody spherical surface (i.e. 4ffi?^), 
the corresponding radius would be half that obtained from a circular region. 



© 2012 RAS, MNRAS OOO-fTiriOl 



Accreting magnetars 5 



Table 4. RGS model parameters. 



Parameter 



ABS*(BB) 



ABS*(BB+PO) 



Nnixlifi^ cm-^) 


0.10±0.02 


0.16±0.04 


kTtt (keV) 


2.1()±0.30 


1.95±0.()3 


Normjfc 


0.0098±0.0027 


0.0082±0.0004 


Xl (dof) 


1.19(702) 


1.19(701) 



Power-law photon index fixed at F = 0.97 

neutron star. This would be consistent with a reprocessing region 
(ring) surrounding the hot spot, the accretion mound or the top of 
the shock front at the impact point and excludes the innermost parts 
of a hypothetical accretion disk. 

Given the limited band pass of XMM-EPIC cameras, one 
might be worried about what fraction of the spectrum is being re- 
ally accounted for by Comptonization. As can be seen in Fig. [3] 
the BMC component is fitting the high energy part of the spec- 
trum while the blackbody is contributing, mainly, at energies be- 
low 3 keV. Nonetheless, it is perhaps convenient to give the fitting 
in terms of a phenomenological model composed by the addition 
of a power law and a blackbody modified at low energies by photo- 
electric absoiption. These parameters are given in the middle part 
of Table [3] As can be seen, the fit is acceptable, albeit with a higher 
value of the reduced ^^. The most noticeable change is the increase 
in temperature of the blackbody component which is now contribut- 
ing at higher energies (peak at 4 keV). The radius of the soft photon 
source would be Rtb — 0.7 km (0.37 km if an spherical emitting re- 
gion is considered), smaller than, but of the order of, that found 
under the Comptonization interpretation. 

The radii obtained from the two-blackbody model are also 
consistent with emission from the polar caps (Table |3}, although 
in this case the model shows the largest residuals. 

The error in /?bb was obtained by propagating the errors from 
the best-fit blackbody normalisation. A', {Rhh = VnDiq)) as given 
by the x^ fit a nd the error on the distance, which was estimated 
to be 0.3 kpc ( IB lav et alj|200q) . Diq is the distance to the source 
normalised to 10 kpc. 

A high signal-to-noise spectrum at low energies is also partic- 
ularly interesting because it allows the determination of the hydro- 
gen column density, Nh- The energy range 0.1-2.0 keV is the most 
suitable band for the derivation of the Nh since cold gas mainly 
absorbs X-ray photons with energies lower than 2-3 keV. We find 
Nh = (0.36±0.01)x 10^^, which is the most accurate determination 
of the absorption to the source measured so far. The Nh allows us to 
derive the opti c al ex tinction to the system. Using the calibration by 
iGiiver & O zel (2009). we obtain Ay = 1.63 ± 0.08 mag or, assum- 
ing a standard reddening law, E(B - V) = Av/3.1 = 0.53 mag, in 
excellent agreement with the value of 0.54±0.05 mag derived from 
optical observations ( IBlav et alj2006r) . This value agrees very well 
with the reddening derived using diffuse interstellar bands, hence 
it is consistent with Galactic absoiption along the line of sight to 
4U 2206+54. We conclude that there is very little matter locally. 



4.2 RGS spectrum 

The low value of the reddening is further supported by the analy- 
sis of the RGS data. In Fig. |4] we show the RGS spectrum of the 
source in the wavelength range 6 to 20 A (= 0.6 to 2 keV). The 
data from RGS 1 and 2 have been added and the resulting spectrum 




12 14 

Wavelength (A) 



Figure 4. RGS spectrum. 



Table 5. Emission line parameters. 



Line 



(keV) 



(keV) 



Flux 

(xlO-^) 

ph s"' cm" 



EW 



(eV) 



FeKor 

Fe XXVI Ly a 

Si XIV Ly a 



6.46± 0.05 
6.99± 0.02 
2.22±0.05 



0.14±0.06 
0.14 
0.05 



8.9±3.0 
5.2±2.0 
6.4±2.0 



25.0±8.4 

16.5±13.2 

7.4±2.5 



Numbers without errors were fixed at the quoted values 

rebinned so that each bin has a minimum signal-to-noise ratio of 7. 
In order to detect the presence of any possible lines, the modeling 
of the underlying continuum is critical. Therefore we take advan- 
tage of the broader energy range of PN and use the parameters de- 
duced in the previous section as input. In this wavelength range the 
main contributor is the blackbody. Consequently we tried a single 
absorbed blackbody as well as an absorbed power law plus black- 
body. In the latter case, the power law parameters were frozen while 
those of the blackbody were let to vary freely. The best-fit parame- 
ters are given in Table |4] Two important conclusions can be drawn 
from this analysis. First, no emission lines from an ionized plasma 
are detected. This is in agreement with the PN analysis where the 
emission line candidates from photoionized species, if present, are 
extremely weak (see Sect. [43}. Second, within the RGS interval 
the main indicators of the photoelectric absorption are the Ne edge 
at 14.2 A and the Fe L edge at 17.5 A. As can be seen, the edges are 
extremely shallow. This supports the low value of the photoelectric 
absorption obtained from the analysis of PN data. 



4.3 Detection of emission lines 

The X-ray spectrum of 4U2206+54 seems to be devoid of any 
strong emission features. The high sensitivity of XMM, however, 
allowed us to detect, for the first time, weak emission lines. The 
residuals are clearly structured around the Fe complex region, at 6.5 
keV and also around 2 keV. The region around 6.5 keV is shown 
in Fig. 6. These features can be fitted with two gaussians centred 
at 6.46±0.05 keV and at 6.99±0.02 keV respectively. Both lines 
show a width of tr = 0.14 keV but only in the first case the er- 
rors could be constrained (±0.06) while we fixed the width for the 
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Figure 5. Residuals in number of sigmas of tlie fit to EPIC-PN energy spec- 
trum around tlie iron complex region. The lower and upper panels show the 
residuals with and without the inclusion of two Gaussian profiles, respec- 
tively 



second one. The lines are significant at three sigma. When the gaus- 
sians are added, the;^f^ of the fit decreases from 251 for 202 degrees 
of freedom to 214 for 198 degrees of freedom. The energy of these 
components is consistent with the Fe Ko- (6.4 keV) transition from 
neutral Fe and Fe XXVI Lya (6.96 keV) transition from ionized 
H-like Fe. The parameters of these lines are shown in Table |5] As 
can be seen , even though the detection is firm, the equivalent width 
is very small (25 eV and 16.5 eV, respectively). Consistently with 
the extreme weakness of the Fe Ka line, no Fe edge at 7.14 keV is 
detected. 

The residuals around the 2 keV region can be fitted with a 
narrow Gaussian at 2.22±0.05 keV. The centroid energy is consis- 
tent with the Si XIV Lya transition of H-like Si. The line appears 
very narrow and can not be resolved by the instrument. We fixed 
the width to 0.05 keV The significance of this line is, however, low 
(2a). 



5 THE ISGRI SPECTRUM 



Previous works dTorreion et all2o64t iMasetti et al.l2004 iBlav et all 



120051 : IWang l2009f) have reported the possible presence of a cy- 
clotron resonant scattering feature (CRSF) at around 30 keV. All 
the detections of this feature were tentative because the statisti- 
cal significance of any individual detection was <2cr. Since the 
gaimna-ray imager on board INTEGRAL is sensitive to X-rays in 
the range 15-1000 keV and has been in operation for nearly 12 
years, it is a promising instrument to search for the cyclotron line. 
We obtained an ISGRI spectrum for each one of the SCW for which 
the source appeared with a detection level above la and lay within 
10 degrees of the centre of the pointing. An average spectrum was 
created from the 37 resulting spectra. The average spectrum ob- 
tained from the 37 SCWs does not show any significant feature at 
the energy of the previous reported line. A single power law com- 
ponent gives an acceptable fit ix' = 28 for 21 degrees of freedom) 
in the range 20-180 keV The best-fit photon index is 2.65 ± 0.08. 
Figure |6] shows the combined EPIC-PN and ISGRI spectra. 
The blackbody plus BMC and two narrow gaussians model gives a 
very good fit (reduced ;^'^ « 1). The best-fit values of the spectral 
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Figure 6. Joined EPIC-PN and ISGRI spectra. The dotted lines correspond 
to the blackbody plus bulk comptonization model of Table |3] 



parameters obtained from the fit of this model to the joined EPIC- 
ISGRI spectrum are consistent with those in Table[3] within errors. 
The energy spectral index is slightly larger, 1.7, compared to 1.5 of 
the EPIC-PN spectrum. 



6 DISCUSSION 

4U 2206+54 is a persistent X-ray source with relatively low lumi- 
nosity (Lx ^ 10'*) compared to other accreting X-ray pulsars. It 
shows X-ray va riability on timescales of a few days with an am- 
plitude of ~ 5 JReig et alJi 2009) and on timescales of years by a 
factor of <100 JRibo et al.li2006 ). The short-term variability has 
most likely an orbital origin (a varying mass accretion rate due to 
its elliptical orbit, e = 0.4), while the long-term variability is prob- 
ably due to osci llations in the wind density from the donor star 
JRibo et alj|2006l) . 

Considerable progress has been made on the nature of this 
enigmatic system in the past years. The optical companion is an 
09.5Vp, where the p indicates that it is a peculiar star. The pecu- 
liarity of 4 U 2206-1-54 is that the helium lines observed in the opti- 
cal sp ectra JBlav et al.ll2006f) and stellar wind velocity dRibo et al.l 
120061) are stronger and slower, respectively, than what is expected 
for a normal main-sequence 09.5 star. The highly debated issue of 
the nature of the compact companion, whether it is a neutron star 
or black hole, has been definitively sol ved in favour of the former 
with the discovery of X-ray pulsations (IReig et al.l2009l) . 

4U 2206-1-54 has never been observed by any high-sensitivity 
instrument below 2 keV (BeppoSAX and Siizaku have ~ 10 times 
less effective area at ~ 0.3keV and 1.5 keV, respectively than 
XMM-Newton). Thus the XMM-Newton observations represent an 
excellent opportunity to investigate the emission properties of 
4U 2206+54 in the soft X-ray band. 



6.1 The source environment 

Fluorescent iron line emission provides strong evidence for the 
presence of material in the vicinity of the X-ray source. All pre- 
vious observations of 4U 2206+54 failed to detect any iron line 
at or near 6.4 keV |Negueruela & ReidlJOOll: iTorreion et ^|2004 



iMasetti et alJl2004lReig et alj |2009i). Upper limits for narrow (0.1 
keV) and broad (0.5 keV) lines had been set in < 56 eV and < 65 
eV, respectively JMasetti et alj2004l) . These results imply that there 
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is not much material intrinsic to the system. The agreement be- 
tween the_coloui- excess E(B - V) derived from optical observa- 
tions JBlav et ani2005h and that obtained from our XMM-Newton 

observations (Sect l4.lt would support this result. 

^In a recent survey of HMXB with Chandra, JTorreion et al.l 

I2OIOI) showed that the fluorescence Fe Ka line is present in virtu- 
ally 100% of the wind-fed HMXBs studied. The curve of growth 
shows that the reprocessing material is a spherical shell of near 
neutral Fe, formed in the stellar wind. It is surprising that although 
4U 2206-1-54 is a wind-fed system this line has never been observed. 
The higher spectral resolution and sensitivity of XMM-Newton be- 
low 10 keV compared to previous missions allowed us to report the 
first detection of emission features at 6.4 keV and 7.0 keV. These 
features can be interpreted as Fe Ka and Fe XXVI Lya from pho- 
toionized plasma, respectively. 

Although statistically significant, these lines are weak. Some- 
how the compact object is not illuminating the stellar wind, as in 
the typical supergiant systems, even though the accretion process is 
taking place, efficiently. The source shows random variations in the 
light curve typical of wind accretors, as expected for a high-mass 
X-ray binary. Our results confirm that the X ray source, however, 
is not deeply embedded into the stellar wind of the mass donor, yet 
it is continuously accreting. The large orbital separation would im- 
ply that the wind is too tenuous when it reaches the magnetosphere 
but the fraction of the stellar-wind mass-loss rate captured by the 
neutron star would still b e large enou gh to power the X-rays. This 
fraction is (see e.g iBozzo et alj2008l) 

where a is the orbital separation and i?acc the accretion radius, that is 
the distance at which the inflowing matter is gravitationally focused 
toward the neutron star. It is given by 

ifacc = 2GMIvl = 3.7 X 10'^ v^^ cm 

The stellar wi nd in 4U 2206-H5 4 is abnormally slow with v = 
3.5 X 10^ cm s"' JRibo et ahluOOfih . Thus the accretion radius in 
4U 2206-H54 is i?,,, = 3 x 10" cm 

The orbital period of 4U 2206-1-54 is still not known. Two 
modulations in the long-term X-ray light c urve at 9.6 d and 19.2 
d have been attributed to the orbital period dCorbet et al.COOTh . In 
either case, the neutron star is orbiting at a very large distance from 
the donor (~ 7R. if />„* = 9.7 d and ~ lOR, if P„a = 19 d). As- 
sumin g the radius of 4U 2206-1-54 to be R, = 7.3 Rg JRibo et al.l 
[200g), the orbital separation would result in a = 3.6 x 10'^ cm and 
a = 5.3 X 10'^ cm for the short and long orbital period solutions, 
respectively. 

Thus the fraction of captured material would be Mcupt/M^y ^ 
0.002 and 0.001, respec tively. The mass- loss rate of an 09.5V star 
is ~ 3 X 10->* Mo yr-i teibo et al.ll200^ . Thus M^apt ^ 3 x lO"" 
Mq yr"' or 2 X 10" g s"', which is enough to power a luminosity 
of ~ La- a GMNsMcapt/J^NS ~ 3 X 10^' erg s"', as observed. 

6.2 The X-ray spectrum 

In order to investigate the nature of 4U 2206-1-54, it is illustra- 
tive to compare its X-ray spectrum with that of other types of 
strongly magnetised neutron-star systems: accreting X-ray pulsars 
and anomalous X-ray pulsars (AXP). The former contain neutron 
stars with magnetic fields B ~ 10'^ - 10'^ G, while AXPs are 
thought to harbour magnetars with B ~ 10'^ - 10" G. 

The spectra of many accreting X-ray pulsars show excess 
emission below ~2 keV. When this soft excess is modeled with 



a blackbody component, the temperature and emis sion radius 
strongly depend on the source luminosity dHickox et al. 20041). 
In Be/X-ray binaries at low luminosity {Lx ^ 10'^ erg s"'), this 
soft component is believed to represent emission from the po- 
lar caps. The main piece of evidence supporting this interpreta- 
tion is the size and temperature of the emittin g area. The radius 
of th e emitting region is typically < 0.5 km I Mukheriee & Paul 



:til2l 



2005:: ILa Palombara & Mereghettil2"006l.l2007l;lLa Palombara et al. 
20091) . significantly smaller than the assumed canonical size of the 
neutron star (~ 10 km). Thus, isotropic emission from the surface 
of the neutron star is ruled out. The blackbody temperature is > 1 
keV. 

The X-ray spectra of AXPs in the 0.5-10 keV band are very 
soft and can best be described by a blackbod y {kT ~ 0.5 keV) plus 
a very steep power -law model with F ~ 3 - 4 jWoods & Thompson! 
120061 : lKaspil2007n or by two blackbody components w ith charac- 
teristic temperatures <1 keV JHalpem & Gotthelj|2005l) . 

The hard X-ray spectra (above 20 keV) of AXPs is charac- 
terised by a hard power law tail extending up to 150 keV without 
any sign of turn over. The power-law photon index range F ~ 1-2, 
that is, much ha rder than the power laws seen below 10 keV 
JGotz et alJI 20061) . The X-ray spectrum of accreting X-ray pulsars 
above 20 keV also exhibits power law tails. However, these are 
significantly softer (F ~ 2 - 3) than those of AXPs. Als o, many ac- 
cretin g pulsars show exponential cutoffs at 30-50 keV JBlav et al.l 
[20T3). 

The spectrum of 4U 2206-1-54 shows a soft excess at low en- 
ergies. Irrespectively of the model used, we find that the size of 
the emitting area in 4U2206H-54 is incompatible with the pres- 
ence of an accretion disc and fits in the polar cap scenario. The 
temperature, kT\^\, = 1.63 keV and radius, R\^\, = 370 m of the 
blackbody component (Table [3} and the photon index (F = 0.94) 
are c omparable to that seen in o t her low-luminosi ty accreting pul- 
sars ( La Palombara et al.l l2009t 'Reig et al. 2009, and references 
therein). However, unlike most accreting X-ray pulsars, but sim- 
ilarly to AXPs, the hard power-law tail in 4U2206-t-54 does not 
show any evidence for a cutoff below ~ 150 keV (Fig.|6l. The hard 
spectrum in 4U 2206-1-54 is not as hard as in AXPs, though. 

6.3 Is 4U 2206-H54 a magnetar? 

4U2206-f54 has been observed by many space observatories. A 
cyclotron resonant scattering featu re (CRSF) at ~30 k eV has been 
reported in observati ons by RXTE JTorreion et all2004l). BeppoSAX 
JMasetti et al.1l2004h and INTEGRAL telav et al.l l20^ . All these 
detections were marginal (at <2cr). Wang (2009) analysed INTE- 
GRAL data of 4U 2206-1-54 and reported a harmonic of this feature 
at ~60 keV, although with very low significance. The alleged CRSF 
was only detected during a high state {Lx ~ 2x 10^' erg s"') but not 
during a low state {Lx « 4 x 10^"* erg s"') in the data analysed by 
Wang (2009). However, just a few days after the high state where 
the detection of the harmonic was reported, 4U 2206-1-54 exhibited 
two short-lived hard flares with peak luminosity Lx » 1.3 x 10^^ 
erg s ' and no CRSF was observed dWanduOlOl) . Quoted values 
of the l uminosity are for the 20-10 keV r ange and a distance of 
2.6 kpc dBlav et alj|2006l) . Likewise. JReige t al. (2009) did not find 
any evidence for a cyclotron line in the RXTE observations that 
led to the discovery of the X-ray pulsations. The observations by 
iReig et al.l (2009) are particularly constraining because they anal- 
ysed more than 138 ks worth of data with instruments capable to 
detect a cyclotron line at 30 keV. The fact that this feature, if real, 
is transient argues against its interpretation as a CRSF. 
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On the other han d, the cyclotron line shape is a strong func- 
tion of viewing angle JArava-Gochez et alJuOOOl) . Because of the 
long pulse period of 4U 2206+54, it might be possible that the 
cyclotron line in 4U 2206+54 is enhanced at certain pulse phases 
and suppressed at others. When analysing a pulse-phase averaged 
spectrum, the significance of the detection would depend on how 
many spectra at the right pulse phase are used. We explored this 
idea by performing pulse-phas e resolved spectrosco py on the 2001 
RXTE observations analysed in lTorreion et all ( 120041) and on ISGRI 
data taken on Rev. 383 (December 2005) and Rev. 510 (December 
2006), but we did not find any significant change with pulse phase. 
This kind of analysis is hampered by low statistics. 

An independent estimate of the magnetic field strength can be 
obtained considering the evolutionary track of the neutron star's 
spin period. Finger et al. (2010) derived an average rate of the spin 
frequency change between a BeppoSAX (November 1998) and a 
Suzakii (May 2007) observation of y = (-1.7 ± 0.3) x lO"''* Hz 
s"' , indicating that the rotation of the neutron star is slowing down. 
Based on the slow pulsation period and the spin down rate, they 
show using a magneto-rotation model for the evolution of neutron 
stars, that the 5500 s pulsations in the X-ray flux of 4U 2206+54 can 
be explained provided the neutron star in this system is a magnetar 
whose surface field at the present epoch exceeds 10'* G. Indeed, 
the spin down rate is related to the magnetic moment of the neutron 
star by (see e.g Ikhsanov & Besk r ovn aya 2 010) 

A<,„ « 1032^7"-(M,.5)'/^/;f vif^Pjsoo G cm^ 

where where 745 = //lO'*' g cm' is the star's moment of iner- 
tia, k, is a dimensionless parameter that accounts for the geome- 
try of the accretion flux and it is close to 1 for spherical accretion, 
v_i4 = v/lO"''' Hz s"' is the rate at which the neutron star rotation 
changes with time, and M is the mass of the neutron star in units 
of 1.5 Mq. From a detailed analysis of th e magneto-rotational evo- 
lution of the neutron star in 4U 2206+54. llkhsanov & Beskrovnaval 
( 1201 oh concluded that its surface magnetic field may range from 
5x 10'^ G to 3x 10'^ G, depending on the geometry of the accretion 
flux outside the magnetosphere, with a most probable value in the 
interval (6 - 9) x 10''' G. This scenario also requires that accretion 
in 4U 2206+54 must occur without the formation of an accretion 
disc. If an accretion disc were present then the observed evolution 
of the rotation period in 4U 2206+54 would imply a magnetic field 
strength larger than 3 x 1 0'^ G, which it is r ather difficult to account 
for with current models (llkhsanov & Beskrovnaya 2010) . The diffi- 
culty to accommoda te an accretion disc was already mentioned by 
lTorrei6net"^ ( l2004 . 

The estimate of the magnet i c field strength based on the spin- 
down rate given by Finger et al.l (120101) was determined using only 
two observations. As the authors state, "it is very unlikely that the 
spin-down rate obtained from the BeppoSAX and Suzaku observa- 
tions is valid on the long term since the pulsar would reach zero 
frequency in 300 yr". Suiprisingly, the spin-down rate obtained us- 
ing all publicly available measurements of the spin period and the 
new value f r om XM M-Newton agrees well with that reported by 
iFinger et al.l ( 120100 . As can be seen in Fig |2] the spin-down rate 
does not show evidence of a break-up. A linear fit to all the mea- 
surements of the spin period of 4U 2206+54 gives a spin-down rate 
of V = (-1.5 ± 0.2) X lO"'"* Hz s"' with a correlation coefficient of 
0.97. This rate implies a magnetic field of the order of 6 ~ 10'* G. 

If the 1998 BeppoSAX data point, i.e., the oldest measurement, 
is removed from the fit, then the spin-down rate decreases to half 
the previous value, v = (-7.5 ± 2.0) x 10"'^ Hz s"', which would 
indicate that neutron star has started to rotate at a faster rate. Still, 



the magnetic field derived from that spin-down rate (see equation 
above) would be more than an order of magnitude stronger than 
common values seen in accreting X-ray pulsars. 

Recently, a new theory of quasi-spherical accr etion in X- 
ray pulsars have been developed (IShakura et al.ll2012l) . where the 
dipole fiel d in a win d-fed neutron star is given by (see e.g. 
iPopov & Turollal2012h 
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where M^, is the mass accretion rate in units of 10'* gr s"', V300 is 
the wind velocity in units of 300 km s"', Piooo is the spin period 
in units of 1000 seconds and Porb.300 is the orbital period in units 
of 300 days. Substituting the values measured for 4U 2206+54: 
Mif, = 0.2, V300 = 1-17, Piooo = 5.5, and P^-i.joo = 0.064, the 
resulting magnetic field is S ~ 1.6 x 10'* G, in agreement with the 
magnetic field obtained from the spin-down rate. 

For a dipole magnetic field of ~ 10'* G, the electron cyclotron 
line would appear at £ > 500 keV, where the photon statistics are 
not good enough. On the other hand, a proton cyclotron line would 
appear at £ ~ 0.5(fi/10'*G) = 0.5 keV. Although a line with this 
energy should be observable with XMM-Newton detectors, it is in 
a region affected by strong interstellar absorption. Furthermore, no 
significant lines have been d etected in the persistent emission of 
magnetars ( lMereghettill2008h . 

The picture that emerges of 4U 2206+54 from the available 
data at various bands of the electromagnetic spectrum is that of an 
intriguing system: a highly magnetized (6 ~ 10'* G) neutron star 
accreting matter from the stellar wind of a 0-type main-sequence 
star. Is 4U 2206+54 a magnetar? 4U 2206+54 does not fit in the tra- 
ditional definition of magnetars, namely, very strongly magnetized 
(B ~ 10'* — 10'^ G), isolated neutron stars powered by magnetic 
energy. First, the magnetic field in 4U 2206+54 is in the lower end 
and most likely below the canonical range of magnetic fields in 
magnetars. Second, it is part of a binary system. Third, X-rays are 
powered by accretion and the contribution to the X-ray luminosity 
from dissipation of the magnetic field is expected to be two to three 
orders of magnitude lower. 

The characteristics of 4U 2206+54 would be closer to, but 
not fully consistent wit h the definition of anomalous X-ray pulsars 
dMereghetti et alj|2002h : a spinning down pulsar, with a soft X-ray 
spectrum, apparently not powered by accretion from a companion 
star, and with a luminosity larger than the available rotational en- 
ergy loss of a neutron star. As shown in Sect. 16.21 the X-ray spec- 
trum of 4U 2206+54 difi'ers from that of AXPs. 

Alternatively, one can simply define a magnetar as a neu- 
tron star with a magnetic field strength larger than the quan- 
tum critical value B, = m^cV3/he = 4.4 x 10'^ G at which 
the energy between Landau levels of electrons equal their rest 
mass. In this case, 4U 2206+54 would contain a magnetar ac- 
creting matter from a luminous companion. A more appropri- 
ate name for 4U 2206+54 wo uld then be accreting magnetar 
( llkhsanov & BeskrovnavaluOlOh . 



6.4 Accreting magnetars as a class of systems 

The idea that wind-fed X-ray pulsars may contain super critical 
highly magnetised neutron star has been discussed in previous 
works. The supergiant X-ray binary 2S 0114+650 is the accret- 
ing X-ray pulsar with the longes t period (2.7 hr). No cyclotro n line 
has been reported in this system. iLi & van den Heuvell ( ll999l) argue 
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that, to interpret such a long spin period, the magnetic field strength 
of th is pulsar must have been initially >10''* G. 

IBozzo et alj ( 120080 showed that the large luminosity swings 
between quiescence and outbursts (up to ~ 10^) seen in super- 
giant fast X-ray transients (SFXT) can be explained as the result of 
transitions across different inhibition of accretion regimes. These 
regimes include inhibition due to the magnetic and centrifugal bar- 
riers. In either case, magnetar-like fields can be attained if the spin 
period of the system is long (>1000 s) or the system is luminous 
(Lx > 10^'^ erg s"'). 

Another candidate for accreting ma gnetar is the wind -fed su- 
pergiant X-ray binary GX 301-2. D oroshenko et al.lteOlOl) showed 
that a very strong magnetic field (2 - 3 x lO''' G) is required to 
explain the long pulse period of this system. A cyclotron line at 
30—45 keV, implying a magnetic field of the ord er 5 x 10'^ G 
JKrevkenbohm et alj |2004| : iLa Barbera et alj |2005|) has been ob- 
served in GX 301-2. This magnetic field is well below the criti- 
cal value. However, these contradicting results can be reconciled if 
the line-forming region is situated far above (2 - 3 j?ns ) the neu- 
tron s tar's surface, i.e., in the accretion column (Doroshe nko et al.l 
|201(J) . Another system that might belong to the class of highly- 
magnetised HMX B is IGR J 1635 8-4726, although its nature is 
much less certain jPatel et al.ll2007r) 

The main difference of 4U 2206-1-54 with respect to 2S 
0114+650, GX 301-2 and SFXTs is that 4U 2206+54 contains a 
main-sequence star, hence a still younger system than those with 
supergiant companions. 4U 2206+54 would be the n atural pro- 
genitor of supergiant X-ray binaries (IRibo et al.ll2006h . However, 
4U 2206+54 may not be the only non-supergiant system with a 
likely highly-magnetised neutron star. The magnetic field of the 
Be/X-ray binary X Per ma y also lie in the magnetar range if a slow 
accretion wind is assumed JDoroshenko et al.ll2012l) . 

By comparing the common properties of these systems we can 
identify the following general characteristics of accreting magnetar 
binaries: 

- They contain slowly rotating neutron stars (Pspin ^ 500 s) with 
magnetic fields above the quantum critical value (4.4 x 10'^ G) but 
lower than a few times 10'* G. 

- They are persistent X-ray sources with luminosities in the 
range 10''^ - 10^* erg s"' with long-term amplitude of X-ray vari- 
ability below two orders of magnitude. 

- They are young systems (<10^ yr), undergoing spherical wind 
accretion from luminous companions (OB stars) with abnormally 
low stellar wind velocity (<400 km s"'). 

- Although brief spin-up episodes may occur, spin down domi- 
nates over the long-term with relatively short spin-down time scales 
T ~ P/P < 10=* yr. 



7 CONCLUSION 

With the discovery of X-ray pulsations, the highly debated issue 
of the nature of the compact companion has been solved, but have 
posed new questions. The combination of slow pulsations and long- 
term spin down of the neutron star can be understood assuming 
that the magnetic field of the neutron star is of the order of lO''' 
G. We have obtained a new value of the spin-down rate and shown 
that it is consistent with a supercritical magnetic field. We searched 
for the presence of a cyclotron line in the hard X-ray spectrum of 
4U 2206+54, but we did not find any significant feature at the pre- 
viously reported energies. 4U 2206+54 adds up to a growing pop- 
ulation of highly-magnetised neutron stars in binary systems, with 



magnetic field strengths in the magnetar range, whose X-ray emis- 
sion is powered by accretion. 
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